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Abstract. We show how an operation of inf-convolution can be used to ap- 
proximate convex functions with C 1 smooth convex functions on Riemannian 
manifolds with nonpositive curvature (in a manner that not only is explicit but 
also preserves some other properties of the original functions, such as ordering, 
symmetries, infima and sets of minimizers), and we give some applications. 



1. Introduction and main results 

Smooth approximation is an old subject. Its importance lies on the fact that 
most analytical tools for studying the properties of functions defined on a normed 
space or on a Riemannian manifold require some degree of differentiability of the 
considered functions. However, many functions which arise naturally from geomet- 
rical problems on manifolds are only continuous (or even merely lower semicontin- 
uous). One is thus tempted to approximate those functions by smooth functions to 
which one can apply more powerful analytical methods and obtain some informa- 
tion about the behavior of the approximations which will hopefully be shared with 
the original (nonsmooth) functions. 

The theory of convex functions is also an old subject which plays an important 
role in many areas of mathematics. In Riemannian geometry it has been used, for 
instance, in the investigation of the structure of noncompact manifolds of positive 
curvature by Cheeger, Gromoll, Greene, Meyer, Siohama, Wu and others, see ^3 
EU EB EH EH IH (recall that a function / : M -> R defined on a Riemannian 
manifold M is said to be convex provided the function R 9 t i— ► / o 'y(i) £ R is 
convex for every geodesic 7 on M). The existence of global convex functions on 
a Riemannian manifold has strong geometrical and topological implications; for 
instance it is shown in |l(Jj that every two-dimensional manifold which admits a 
global convex function which is locally nonconstant must be diffeomorphic to the 
plane, the cylinder, or the open Mobius strip. 

Along with the papers cited above and the references therein, we must mention 
the important work of Bangert's on convex sets and convex functions on Riemannian 
manifolds, see 0EJE]; he showed in particular that Alexandroff 's Theorem about 
almost everywhere twice differentiability of convex functions on R n can be extended 
to convex functions on finite-dimensional Riemannian manifolds (providing as well 
a smart proof of Alexandroff 's theorem on R"). 
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The aim of the present paper is to study to what extent one of the most useful 
methods for regularizing convex functions on normed spaces, namely that of infimal 
convolution, can be successfully used in the setting of Riemannian manifolds. 

Let us first have a quick look at the three main methods (that is, partitions of 
unity, integral convolution with a sequence of mollifiers, and inf and sup convolution 
formulae) that are used in normed spaces to approximate continuous functions by 
smooth functions, and see how they can be adapted to the case when one wants to 
regularize a convex function / defined on a Riemannian manifold M. 

Partitions of unity are useless in this setting because, even in the case when 
M = R", they do not preserve convexity of the function /. 

The integral convolution with a sequence of mollifiers reveals itself as the perfect 
tool when M = R™, because the integral convolution of a convex function / with 
any integrable function g with compact support, that is, 



is a convex function. In El EU Greene and Wu studied to what extent those 
integral convolutions with mollifiers can be used to regularize convex functions 
defined on finite-dimensional Riemannian manifolds M (and applied this method 
to prove several theorems about the structure of complete noncompact manifolds 
of positive curvature). It turns out that this method works out in Riemannian 
manifolds only when the original function / is strictly convex. More precisely, let 
k : R — > R be a nonnegative C°° function with support in [—1,1], constant on a 
neighborhood of and satisfying J R „ k(||x||)<Az; = 1, and let us define the functions 



where dfi p is the measure on the tangent space TM p obtained from the Riemannian 
metric of M. Greene and Wu showed that if / : M — > R is a convex function defined 
on an rt-dimensional Riemannian manifold M and K is a compact subset of M, 
then there exists an open neighborhood of K and an £o > such that the functions 
(p e : U — ► R defined above are C°° smooth, converge to / uniformly on K as e — ► 0, 
and are approximately convex in the sense that 



where the infimum is taken over all geodesies j(t) parameterized by arc length and 
with 7 (0) e K. 

Now, a C°° function cp : M — ► R is called strictly convex if its second derivative 
along any geodesic is strictly positive everywhere on the geodesic. A (not necessarily 
smooth) function / : M — > R is then said to be strictly convex provided that for 
every p € M and every C°° strictly convex function tp defined on a neighborhood 
of p there is some e > such that / — eip is convex on the neighborhood. With 
this terminology, the above result implies that if / : M — > R is strictly convex then 
for every compact subset K of M there exists a sequence of strictly convex C°° 
functions (p n = tp Sn such that / = linin^oo (p n uniformly on K. 

However, as Greene and Wu pointed out, this method cannot be used when / 
is not strictly convex, and the problem of approximating (not necessarily strictly) 
convex functions by smooth convex functions on Riemannian manifolds is open. 





liminf(inf^ £ ( 7W )| t J>0, 
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That is one main limitation of the integral convolution technique on Riemannian 
manifolds. Another drawback of this method is the fact that it does not apply to 
functions defined on infinite-dimensional manifolds (even in the case when M is the 
Hilbert space). 

We are left with the third method: infimal convolution. It is well known that if 
/ : X 1U {+00} is a lower semicontinuous convex function defined in X = W l 
or in any infinite-dimensional reflexive Banach space X (such as, for instance, the 
separable Hilbert space), then the inf-convolution formula 

f k (x) := inf{/(y) + fc||z-u|| 2 }, 

yex 

where || • || is any equivalent differentiable norm in X, defines a sequence of C 1 
smooth convex functions fu which converge to / as k — > +00 (uniformly on bounded 
sets if / is uniformly continuous on bounded sets). In fact, a clever combination 
of inf and sup convolutions allows to show that if / is a (not necessarily convex) 
function which is uniformly continuous on bounded sets of a superreflexive Banach 
space then / can be approximated by C 1 smooth functions with uniformly continu- 
ous derivatives uniformly on bounded sets; this was shown first by Lasry and Lions 
[T5| in the case when X is the Hilbert space and then by Cepedello-Boiso |S] for 
any superreflexive Banach space X. 

In this situation it is natural to ask whether infimal convolution formulae can 
be used to regularize convex functions defined on Riemannian manifolds (cither 
finite or infinite-dimensional). That is the question we try to address in this paper. 
Let us describe the main results that we will show in the following sections. If M 
is a complete Riemannian manifold and d is the geodesic distance on M, for any 
function / : M -^RU {+00}, and for A > we define the function 

f x (x)= M{f(y) + ±-d(x,y) 2 } 

y£M IX 

for every x £ M. In Section 2 we collect some general properties of f\ that do 
not depend on the geometry of M and that we will need to use later in our proofs. 
We show for instance how the inf defining f\(x) can be restricted to a suitable ball 
B(x,r x ), and then use the estimates on the radius r x to show that lim A ^ + f\(x) = 
f(x) pointwise whenever f(x) < +00, and that if / is uniformly continuous on 
bounded sets then fx converges to / uniformly on bounded sets. We also see that 
fx has the same inf and the same set of minimizers as / does, and that fx has 
the same symmetry properties as / (that is, if T : M — » M is an isometry and 
f(Tz) = f{z) for all z, then f x (Tz) = f x {z) for all z). 

In Section 3 we assume that / : M -> K U {+00} is a convex function and we 
study under what conditions on M the functions f\ are convex and C 1 smooth. It 
turns out that some assumptions on the geometry of M are necessary in order that 
the fx be convex and C 1 smooth (see Example 1221 below 1 ) ; in particular we must 
require that the distance function d : M x M — > R, (x,y) <— > d(x,y), be uniformly 
locally convex on bounded sets near the diagonal (see Definition below) . Under 
this assumption we prove that the functions f\ are convex and C 1 smooth on any 
given bounded subset B of M, for all A small enough. Moreover, if the distance d 
is convex on all of M x M then the f\ are convex and C 1 smooth on all of M for 
all A > 0. 

In Section 4 we study the question as to which manifolds satisfy the above 
technical assumption that d is uniformly locally convex on bounded sets near the 
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diagonal (resp. convex on M x M). First, we show that for every Riemannian 
manifold with nonpositive sectional curvature, and with the property that the con- 
vexity radius function of M is strictly positive on bounded sets (such is the case, 
for instance, of all complete finite-dimensional Riemannian manifolds), the distance 
d : M x M — > R is uniformly locally convex on bounded sets near the diagonal. 
Secondly, we note that, for every Cartan-Hadamard manifold (that is, every simply 
connected complete Riemannian manifold of nonpositive sectional curvature), the 
distance d is convex on all of M x M . By combining these facts with the results of 
Sections 2 and 3 we obtain, in the finite-dimensional case (see Corollarvll7lbelow % ): 
If M is a complete finite- dimensional Riemannian manifold with sectional cur- 
vature K < and f : M — ► R is a convex function, then, for every bounded open 
convex set U with compact closure U , there exists A > such that the functions 
fx'-M—yM. are convex and C 1 smooth on U for all A € (0, Xq). Moreover, 

(1) fx converges to f uniformly on U . 



(3) f\ has the same inf and the same set of minimizers as f . 

(4) fx has the same symmetries as f ( that is, if f is invariant with respect to 
an isometry T : M — > M, then so is fx)- 

And, in the case of a Cartan-Hadamard manifold (either finite-dimensional or 
infinite-dimensional, see Corollary II 81 below) : 

If M is a Cartan-Hadamard manifold and f : M — > R U {+00} is a lower- 
semicontinuous convex function, then the functions fx-M^M. are convex and C 1 
smooth on all of M for all A > 0. Moreover, 

(1) If f is uniformly continuous on bounded sets then fx converges to f uni- 
formly on bounded sets. 



(3) fx has the same inf and the same set of minimizers as f . 

(4) fx has the same symmetries as f. 

Finally, in Section 5 we consider some corollaries and applications of the above 
results. We show that if C is a closed convex subset of a Cartan-Hadamard manifold 
then the distance function to C, x 1— ► d(x,C) — inf{d(x, y) : y € C}, is C 1 smooth 
on M \ C and the function x 1— > d(x, C) 2 is C 1 smooth and convex on all of M. 
We also note that this result is not true for Riemannian manifolds of positive 
curvature such as the 2-sphere, and therefore the results of Section 3 cannot be 
extended to manifolds of positive curvature. Another consequence is that every 
closed convex subset of a Cartan-Hadamard manifold can be approximated by C 1 
smooth convex bodies of M. Lastly, we note that if M is a Cartan-Hadamard 
manifold and / : M — ► R U {+00} is convex and lower-semicontinuous, then the 
function 



is the unique viscosity solution to the following Hamilton- Jacobi partial differential 



(2) fx<f for all A > 0. 



(2) fx < f for all A > 0. 



u{t,x) := inf {f(y) + -d(x,y) 2 } for t > 0, u(0,x) = f(x) 



equation: 




where u : [0, 00) x M — » R. 
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2. General properties 
Throughout the paper, for a function / :M^KU {+00}, we define 

f x (x)= inf{f(y) + -Ld(x,y) 2 }. 

yeM ZA 

The following Proposition shows how, under certain conditions, the inf defining f\ 
can be localized on a neighborhood of the point x. 

Proposition 1. [Localization] Let M be a Riemannian manifold, f : M — > K U 
{+00} a function satisfying that f(x) > — 1(1 + d(x, x a ) 2 ) for some c > 0, x G M. 
Let x e M be such that f{x) < +00. Then, for all A e (0, ^) and for all p > p, 
where 

2f(x) + c(2d(x,x ) 2 + 1 2 



p = p(x,X,c) : 
we have that 



1 - 2Ac 
1 



f x (x)= inf {f{y) + ^d(x,yf}. 

y€B(x,p) ZA 

Proof. Since 

d(y, x ) 2 < (d(y, x) + d(x, x )f < 2 (d(x, y) 2 + d(x, x Q ) 2 ) , 
we have that 

c (d(x, y) 2 + d(x, x ) 2 ) > -d(y, x ) 2 , 

hence 



that is 



£ (1 + d(y, x ) 2 ) + -^(z, y ) 2 > ( -L _ c ) d ( :B) y )2 _ cd ^ XQ f _ (!) 



2 K Vtf ' uy 7 2A V 2A / v ' a/ v ' u/ 2 

Now, for any given r] > 0, if we set 

r = r(*, A, c, „) := + 2. + c^.xp)' + 1) ^ 1/2 

by using (1) we obtain that, for every y S M with x) > r, 

/(I/) + 7^d{x, y) 2 >-^(l + d(y, x ) 2 ) + ^d(x, y) 2 > 
^-cj d(x, y) 2 - cd(x, x ) 2 - | > 

i - c) ^■Hl.ty^l ) _ „ i(XiX0)2 _ £ _ 

f{ X )+T]> f X (x), 

which implies that 

inf {/(id + ^dfoy) 2 } >/*(*), 

d{y,x)>r LA 

hence 

/ A (x)= inf {/ (j / ) + -ld( a; , 2/ ) 2 }. 

a(y,x)<r LA 
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Finally, since 



lim r(x, A, c, rj) = p(x,X,c), 
77— >o 



it is clear that for every p > p we can find 77 > small enough so that 

P = P(x, A, c) > r(x, A, c, rf) > p(x, A, c), 
and therefore, from the above argument we deduce that 

f x (x)= inf J-d^y) 2 }. 

d{y,x)<p ZA 



□ 



Next we state several interesting properties of this inf-convolution operation, 
such as preservation of order and symmetry properties of the original function. We 
put off studying the conditions under which convexity is preserved until the next 
section. 

Proposition 2. Let M be a Riemannian manifold, f : M — > MU{+oo} a function. 
We have that: 

(1) fx < f for all A > 0. 

(2) If0< X 1 < A 2 thenf X2 < f Xl . 

(3) inf fx = inf / and, moreover, if f is lower semicontinuous then every min- 
imizer of fx is a minimizer of f , and conversely. 

(4) If T is an isometry of M onto M, and f is invariant under T (that is, 
f(Tz) — f(z) for all z € M), then fx is also invariant under T, for all 
A > 0. 

Proof. (1) and (2) are obvious. 

(3) Note that 

inf f x (x) = inf inf {f(y) + ±d(x,yf} = 
¥m ^l{f{y) + ^d{x,yf}= inf f(y). 

(4) We have that 

f x (Tx) = inf {f(y) + ±rd{Tx, y) 2 } = inf {f(Ty) + ^;d(Tx, Tyf} = 

y&M ZA y£M ZA 

inf < {f(Ty) + ±d(x, y) 2 } = m£{f(v) + ^d(x, y) 2 } = f x (x). 

yeM ZA ySJ" ZA 



a 



Now we apply Proposition^to show that, under natural continuity assumptions 
on /, the regularizations fx converge to the original function / as A goes to 0. 

Proposition 3. [Convergence] Let M be a Riemannian manifold, f : M — > M U 
{+00} a function satisfying that f(x) > — 1(1 + d{x, Xq) 2 ) for some c > 0, xq € M, 
and consider 

1 - ^2 ^ 



fx(x)= ■ m l{f(y) + —d(x, y y} 

y£M ZA 



for < A < l/2c. 
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(1) Assume that f is uniformly continuous on bounded subsets of M . Then 
lirnA— >o fx = f uniformly on each bounded subset of M. 

(2) Assume that f is continuous on M . Then lim^o f\ — f uniformly on 
compact subsets of M . 

(3) Assume that f is uniformly continuous and bounded on all of M . Then 
liniA^o fx = f uniformly on M. 

(4) In general ( that is, with no continuity assumptions on f) we have that 
liniA^o f\( x ) — f( x ) f or every x £ M with f(x) < +oo. 

Proof. (1) According to Proposition^ for every x £ M, A <E (0, l/2c), p > p(x, A, c), 
we have that 

f x (x)= inf {f( y ) + ±d(x,y) 2 }. 

a(x,y)<p LA 

Fix R > 0. As is easily shown, a uniformly continuous function on a Riemannian 
manifold is bounded on bounded sets, hence we can find k > so that |/(x)| < k 
for all x £ B{x ,2R). 

For any given e > 0, by uniform continuity of /, there exists 5 > such that if 
y, x £ B(xo,2R) and d(x,y) < S then \ f(x) — f(y)\ < e/3. We can assume 5 < R. 
Now, since 

lim / a + ^ + 

a^o+ V 1 - 2cA / 
there is X e > such that if < A < A e then 



2k + c(2R 2 + l)^ 1/2 



0<P(*.VO<(A <5 
for all x £ B(xq, R), and therefore 

f x (x)= inf {f(y) + ± d (y,x) 2 } 

d(y,x)<5 ZA 

for all x £ B(xq, R), A £ (0, A e ). But, since fx < f for all A, this really means that 
f x (x)= inf {f(y) + -Ld(y,x) 2 }, 

y£A x ZA 



where 



A x := {y £ B(x, 5) : f(y) + ±-d(y, xf < f(x)}. 



By the definition of inf, we can take y x £ A x such that 

/a(x) + | >f(y x ) + ^d(y x ,x) 2 . 

Then, bearing in mind that y x £ A x C B(x, 8) C B(xo, 2R) when x £ B(xo, R), we 
get 

|/(x) - f x {x)\ = f(x) - f x {x) < f(x) - f(y x ) ~ ^d(y x ,x) 2 + | < 

l/W - f(y*)\ + \m f(y x )\ + |^| + | + |= £ 

for all x £ B(xq,R), A £ (0, A e ). This shows that lim /)i ^Q+ fx — f uniformly on 
B(x ,R). 

(2) Let if be a compact subset of M. By compactness, it is easily seen that for 
every e > there exists S > such that if x G K, y G M, and d(x,y) < S then 
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\f{x) — f(y)\ < e/3. One can now repeat the above argument with the precaution 
of always taking x € K c B(xo, R), y € M, y) < 5. 

(3) Choose fc > such that \f(x)\ < k for all x <G M. Let us first observe that the 
inf defining fx can be restricted to the set {y € M : ir) < 2\/fcA} = B(x, 2\/kX). 
Indeed, if d(y, x) > 2\/k\ then 

f(y) + ^d(x, yf > -k + 2k = k> f(x) > f x (x). 

Next, for any given e > 0, the uniform continuity of / provides us with 6 > 
so that \f(y) — f(x)\ < e/3 whenever d(x,y) < S. On the other hand, since 
rim - \_ > o+ 2\/fcA = 0, there exists A e > such that < 2\/fcA < 5 for < A < A £ . 
Then, for any x € M, A e (0, A e ) the inf defining f\(x) can be restricted to the set 
A x ■= {y € B(x,S) : f(y) + ^d(x,y) 2 < f(x)}. Now, the same estimations as in 
(1) above show that \f\(x) — f(x)\ < e. 

(4) is very easy. □ 



3. REGULARIZATION OF CONVEX FUNCTIONS 

In order to see that the operation / — ► f\ preserves convexity we will need to 
use the following Lemmas. 

Lemma 4. Let M be a Riemannian manifold, and F : M x M — > K U {+00} a 
convex function (where M x M is endowed with its natural product Riemannian 
metric). Assume either that M has the property that every two points can be con- 
nected by a geodesic in M, or else that F is continuous and M is complete. Then, 
the function ip : M — > M defined by 

ip(x) = inf F(x,y) 

is also convex. 

Proof. Let 7 : I — > M be a geodesic. We have to see that the function t V(tW) 
is convex, that is, for any to, h € /, and for any s e [0, 1], ip(-f(sti + (1 — s)t )) < 
s^(7(ti)) + (1 — s)'0(7(to)). We may assume (up to an affine change of parameters) 
that to = and t\ = 1, so we have to show that ip(l(t)) < ti>(xi) + (1 — t)i>(xo), 
where x := 7(0) and X\ := 7(1). 

Fix an arbitrary t e [0, 1]. For any e > 0, by the definition of ip, we can pick 
2/o , yi € M so that 

F(xi,yi) < V(#i) + e, and F(x , yo) < ip{xo) + £■ (*) 

Let a : J — > M be a geodesic connecting yo and yi (if such a does not exist then 
we can use continuity of F and Ekeland's approximate Hopf-Rinow type theorem 
to get points yo and y~i close enough to yo and y\ so that (*) remains true when 
replacing yj with y~j , and a geodesic cr connecting yo to yt ; the rest of the argument 
applies without changes). We can also assume that J = [0, 1], yo = °"(0), yi = o"(l). 

It is clear that, because 7 and a are geodesies in M, the path t (■y(t), cr(t)) is 
a geodesic joining the points (xo, yo) and {x\, y\) in the product manifold M x M. 
Now, since 1 1— > F(^(t),a(t)) is convex, we have that 

V(7(*)) - inf F(j(t),y) < F( 7 (t),a(t)) < tF( Xl ,Vi) + (1 - tj^clto) < 

2/GM 

i(V>(zi) +£) + (!- t)(^(a;o) + e) = + (1 - *M*o) + e 
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and this holds for every e > 0, hence we can conclude that t/KtW) ^ to/>(xi) + (1 — 
s)ip(xo). □ 

Lemma 5. Let M be a Riemannian manifold with the property that any two points 
of M can be joined by a minimizing geodesic, and let f : M — > R U {+00} be a 
convex function. Then, for every Xq G M there exists a number c > such that 
fix) > -f (l + d(x, x ) 2 ) for all x G M. 

Proof. Choose ( G D~ f(xo). For a given x S M, let 7 be a minimizing geodesic 
connecting xq to x, say j(t) — exp xa (tv), t 6 [0, d(x, Xo)], for some u € TM Xo with 
llullso = 1- Since / is convex we have f(exp xo (tv)) - f(x Q ) > (G,tv) Xo for every 
t G [0, d(xo, x)], and in particular 

f(x) - f(x ) > (C,d(x ,x)v) Xa > -\\C\\x o d{x ,x)\\v\\ Xo = 
-\\C\\ Xo d(x ,x) > -\\C\\ X0 (l + d(x a ,x) 2 ), 

and therefore 

f(x) > f(x ) - HCllxo (1 + d(x , x) 2 ) > -| (1 + d(x , x) 2 ) 

for all x G M if we put c = 2 (||C||x + l/(»o)|)< □ 

It will be also useful to recall that every convex function / : M — > R which is 
locally bounded is continuous (in fact locally Lipschitz); a proof of this statement 
can be found in |21 Proposition 5.2]. 

In order that f\ is convex whenever / is, we will have to require that the distance 
function d : M x M — ► R is convex on a band around the diagonal of M x M. More 
precisely, we will use the following. 

Definition 6. Let M be a Riemannian manifold. We say that the distance function 
d : M X M — ► R is uniformly locally convex on bounded sets near the diagonal if, 
for every bounded subset B of M there exists r > such that <i is convex on 
B(x,r) x B(x,r), and the set B(x,r) is convex in M, for all x G -B). 

Examples of manifolds satisfying this definition are the cylinder x 2 + y 2 = 1 
in M 3 , the Poincare half-plane, or the subsets of M 3 defined by z = \j(x 2 + y 2 ) 
or z — xy. In general, every complete finite-dimensional Riemannian manifold of 
nonpositivc sectional curvature meets this condition, as we will show in the next 
section. 

Proposition 7. Let M be a Riemannian manifold with the property that any two 
points of M can be joined by a geodesic, and let f : M — ► R U {+00} be a lower- 
semicontinuous convex function. 

(1) Assume that f is bounded on bounded sets and that the distance function 
d : M x M — > R is uniformly locally convex on bounded sets near the 
diagonal. Then, for every bounded subset B of M there exists Xq > such 
that f\ is convex on B for all A G (0, Ao). 

(2) Assume that the distance function d : M x M — > R is convex on all of 
M x M . Then f\ is convex on M for every A > 0. 

Proof. (1) We may well assume B = B(xo,R) for some xq G M, R > 0. Let r > 
be small enough so that the function (x,y) 1— > d(y, x) is convex on B(x, 2r) xB(x, 2r) 
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and B(x, 2r) is convex for every x £ B(xq, R). Let A; be a bound for / on B(xq, 2R). 
We may assume 2r < R. We have that, for every z £ B(xo, 2R), 

hence we can choose Ao > small enough so that p(z, X ,k) < r for all z £ B(xo, 2R) 
and therefore, according to Proposition^ we have that, for every A £ (0, Ao), 

h(z)= inf {f(y) + -L d (z,y) 2 }= inf {/(„) + i-d(*,v) a } 

yeB(z.r) 2A y£B{x,2r) 2X 

whenever z 6 B(x,r), x £ B(xq,R). Note that B(x,2r), as a convex subset of 
M, still has the property that any two of its points can be joined by a geodesic 
inside B(x,2r). Now, assuming < A < Ao, and fixing x £ B(xo,R), because the 
function F{z 1 y) := f(y) + j^d(z,y) 2 is jointly convex on B(x,2r) x B(x,2r), we 
deduce from Lemma0]that z i— » f\(z) — uif y eB(x,2r)F( z ,y) is convex on B(x,r), 
for all < A < Ao- Since x £ B(xq,R) is arbitrary this implies that f\ is convex 
on B(x ,R), for every A £ (0, Ao). 

(2) Here we can use Lemma 0] on all of M x M with no need to localize the inf 
defining f\(x), so it follows that f\ is convex for all A > 0. □ 

Remark 8. Note that in Case (2) of the above proposition we do not require 
continuity of the function /, so / is permitted to take the value +oo at some 
points; in particular we are allowed to take / to be the indicator function of a 
closed convex subset C of M, that is 



5c(x) = 



if x £ C, 
+oo otherwise 



Remark 9. An examination of the above proof and the statement of Lemma 0] 
reveals that, if one assumes that / is continuous and M is complete, it is not 
necessary to require that every two points of M can be connected by a minimizing 
geodesic in M. 

The following Proposition shows that the functions f\ are superdifferentiable at 
a point x if d is differentiablc on a suitable ball around x. We refer the reader to 
for the properties of viscosity subdifferentials on Riemannian manifolds; here we will 
only make use of the very definition of the subdifferential and the superdifferential 
sets of / : M — > [ — oo, oo], namely, 

D~f(x) = {dip(x) : ip £ C 1 (Af, R), / — ip attains a local minimum at x}, 

and 

D + f(x) = {dtp(x) : £ C 1 (M, R), f — ip attains a local maximum at x}, 

the fact that / is differentiable at x if and only if D~ f(x) ^ ^ D + f(x) (in which 
case {df(x)} = D~f(x) = D + f(x)), and that a convex function / : M — ► R is 
everywhere subdifferentiable Theorem 5.3] 

Proposition 10. Suppose that the inf defining f\{x) can be restricted to a ball 
B x = B(x,r x ) of radius r x small enough so that the function y i— * d(y,x) 2 is 
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differentiable on B x , and that this inf is attained at a point y x G B x . Then fx is 
superdifferentiable at x, and 



1 d 

-?d(x,y x )—d{x,y x ) G D + f x {x). 
A ox 



Proof. We have that 



f\{z) ~ fx(x) < f{y x ) + -^d(z 7 y x ) 2 - f(y x ) - -^-d(x,y x ) 2 = 
2^d(z,y x ) 2 - -^rd(x,y x f, 

so 

fx(z) - ^d(z,y x ) 2 < f x (x) - — d(x,y x ) 2 

for every z G B x , that is, fx — jx~d(-,y x ) 2 attains a local maximum at x, hence fx 
is superdifferentiable at x, with d (jx~d(-,y x ) 2 ^) (x) G D + fx{x). □ 

Next we show that convex differentiable functions are automatically of class C 1 . 
In this proof we will make use of the parallel transport of vectors along geodesies. 
Recall that, for a given curve 7 : I — > M, numbers to,ii G /, and a vector Vq G 
TM 7 ( to ), there exists a unique parallel vector field V(t) along -f(t) such that V(to) = 
Vb- Moreover, the mapping defined by Vq >— > V(t\) is a linear isometry between 
the tangent spaces TM 7 ( to ) and TM 7 ( tl ), for each t\ G J. In the case when 7 is 
a minimizing geodesic and 7(^0) — x , 7(^1) — 2/; w e will denote this mapping by 
L xy , and we call it the parallel transport from TM X to TM y along the geodesic 
7. See ^HI for general reference on these topics. The parallel transport allows 
us to measure the length of the "difference" between vectors (or forms) which are 
in different tangent spaces (or in duals of tangent spaces, that is, fibers of the 
cotangent bundle), and do so in a natural way. Indeed, let 7 be a minimizing 
geodesic connecting two points x,y G M, say 7(^0) = £,7(^1) = V- Take vectors 
v G TM X , w G TM y . Then we can define the distance between v and w as the 
number 

||W - Ly X (w)\\ X = \\W - L X y{v)\\y 

(this equality holds because L xy is a linear isometry between the two tangent spaces, 
with inverse L yx ). Since the spaces T*M X and TM X are isometrically identified by 
the formula v — (v, •), we can obviously use the same method to measure distances 
between forms C G T*M X and 77 G T*M y lying on different fibers of the cotangent 
bundle. 

It is also well known that the mapping y 1— > L xy is well defined and continuous on 
a neighborhood of each x € M, in the following sense: if (x n ) converges to x in M 
then expj, {L XXn (v)) converges to exp x (v) uniformly on the set {v G TM X : \\v\\ x < 
6} for some S > (a fact which we use at the end of the proof of the following 
lemma). 

Lemma 11. Let M be a Riemannian manifold, and let f : M — > K fee a differen- 
tiable convex function. Then f is of class C 1 on M . 

Proof. Assume that / is not C 1 , then there are e > 0, a point x G M and a sequence 
(x n ) C M converging to x such that 

\\L XnX [df{x n )]-df(x)\\ x >2s 
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for all n E N. Therefore, for every n E N we can pick h n E TM X with ||/i n |U = 1 
such that 

(L XnX [df{x n )] - df(x), h n ) x > 2e for all neN. 

Since / is diffcrentiable at x, there exists 6 > so that 

/(exp x (tw)) - f(x) - (df(x),tv) x < et 

for all v E TM X with \\v x \\ — 1 and \t\ < S. On the other hand, by convexity of /, 
we have 

(df(x n ),tw) Xn < f(exp Xn (tw)) - f(x n ) 

for all w E TM Xn with ||ttf|| Xn = 1 and \t\ < S. By combining these inequalities we 
get 

2eS < (L XnX [df(x n )] - df(x),h n ) x 5 = (L XnX [df(x n )],Sh n ) x - (df(x),Sh n ) x = 
(df(x n ),SL XXn h n ) Xn - (df(x),5h n ) x < 
f( ex Px n ( SL xx n h n )) - f( x n) + f(x) - f(exp x (5h n )) + ed = 
f(exp Xn (SL XXn h n )) - f(exp x (5h n )) + f(x) - f(x n ) + eS -> + + sS 

(by continuity of /, exp and the parallel translation L), so we get that 2eS < sS, 
which is not possible. □ 

Now we can prove that, under the same assumptions as in Proposition if / is 
convex then f\ is of class C 1 for A > small enough. 

Theorem 12. Let M be a Riemannian manifold and let f : M — > M U {+00} be 
a lower semicontinuous and convex function. Assume that every two points of M 
can be connected by a minimizing geodesic in M . 

(1) Suppose that f is bounded on bounded sets and that the distance function 
d : M x M — > K is uniformly locally convex on bounded sets near the 
diagonal. Then, for every bounded open convex subset B of M there exists 
\q > such that f\ is a C 1 smooth convex function on B, for all A E (0, Ao). 

(2) Suppose that the distance function d : M x M — > R is convex on all of 
M x M . Then f\ is a C 1 smooth convex function on M , for every A > 0. 

Proof. (1) We can give an almost self-contained proof of this in the finite-dimensional 
case, so let us first assume that dimM < +00. We may also assume B = B{xq,R) 
for some xq E M, R > 0. Since the index of injectivity x <— > i(x) is a continuous 
positive function, it is bounded below by a positive number on the compact sub- 
set B(xo,R) of M. This implies that there exists ro > such that the function 
y 1— > d(x,y) 2 is C 1 smooth on B(x,ro) for every x E B(xo,R). We can obviously 
assume that r < r and repeat the argument of the proof of Proposition to get a 
A > such that fx is convex on B(xq, R) for all A E (0, Ao) and, moreover, that, 
for every x E B(xo,R), the inf defining f\{x) can be restricted to the ball B{x,r), 
which is contained in B(x, ro) (so that, in particular, y 1— > d(y, x) 2 is C 1 smooth on 
B(x, r)). Besides, this inf is attained, because the involved functions are continuous 
and the ball B(x, r) is compact. According to Proposition 1 1(JI we then get that fx 
is superdifferentiable at x, for every x E B(xq, R), A E (0, Ao). 

On the other hand, since fx is convex on B(xo, i?), we know that fx is subdiffer- 
entiable on B(xo, R) (see Theorem 5.3]). That is, fx is both subdifferentiable and 
superdifferentiable at each point of B(xo, R), hence fx is differentiable on B(xq, R) 
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(see Proposition 4.6]). Since fx is convex, Lemma ITD allows to conclude that fx 
is C 1 smooth on B(xo, R) for each A € (0, Ao). 

Let us now consider the case when dimM = +00. Since M has the property 
that any two of its points can be connected by a minimizing geodesic in M, and 
fx, being convex on B for all A € (0, Ao), satisfies D~ f\(x) 7^ for all x 6 B, we 
can apply Theorem 11 of to get that fx is differentiable at every point x € B, 
hence (by Lemma fTTjl of class C 1 on B. 

(2) As in case (1), let us first give a self-contained proof for the finite-dimensional 
case. If the distance function d : M X M — > K is convex on all of M x M then 
y 1— ► <i(y, x) is convex on M for all x € M, which implies that the cut locus of x is 
empty for every x € M and that the function y 1— > d(y, x) 2 is differentiable on all 
of M for every x S M. On the other hand, we claim that the inf defining fx(x) is 
attained for every x € M. Indeed, fix x G M with f(x) < +00. From the proof 
of Lemma 13 we know that there exists c = c x > such that f(y) > —cd(x,y) for 
every y £ M. Then we have 

if y) — > +00, so there exists i? > large enough so that if <i(x, y) > R then 

/(y) +!^y) 2 >/(*)> /a 

hence 

A(aO= mf {f(y) + l-d(x,y) 2 }, 

y£B(x,R) 

and now it is clear that this inf is attained because f + -^d(-,x) 2 is lower semicon- 
tinuous and B{x,R) is compact. 

Therefore, according to Proposition ^3 fx is superdifferentiable. Because fx is 
convex, this means that fx is diferentiable, hence C 1 smooth on M, for all A > 0. 

In the infinite-dimensional case we only have to bear in mind that, according 
to Proposition |7f 2) , fx is convex on all of M, so the same proof as in case (1) 
applies. □ 

Remark 13. If one assumes that / is continuous and M is complete, it is not 
necessary to require that every two points of M can be connected by a minimizing 
geodesic in M. 



4. Which manifolds do the above results apply to? 

Let us say a few words about the manifolds satisfying the assumptions of Theo- 
rem^] The following Theorem is a restatement of ^| Theorem IX. 4.3, p. 257] 

Theorem 14. Let M be a Riemannian manifold with seminegative sectional curva- 
ture K < 0, and U a convex open set. Let [3\, P2 be disjoint geodesies in U , defined 
on the same interval. Let a t : [a, b] — * U be the unigue geodesic joining (3\(t) with 
/3 2 (t), and let £(t) = length(a t ), that is, l(t) = d(0i{t), fo(t)). Then I" (t) > for 
all t, and in particular £(t) is a convex function. 
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From this Theorem it is immediate to deduce that the above results on regular- 
ization of convex functions apply to manifolds of seminegative sectional curvature, 
as we next see. 

Corollary 15. Let M be a Riemannian manifold with sectional curvature K < 0. 

(1) Suppose that M has a convexity radius function which is strictly positive on 
bounded subsets of M ( such is the case, for instance, of a complete finite- 
dimensional Riemannian manifold M). Then the distance function d is 
uniformly locally convex on bounded sets near the diagonal of M x M. 

(2) Suppose that M is simply connected (which, together with the curvature 
assumption, amounts to saying that M is a Cartan-Hadamard manifold). 
Then the distance function d is convex on all of M . 

Proof. (1) Let B be a compact subset of M. Since the convexity radius function 
x i ► c{x) is bounded below on B by a number r > 0, we have that the open ball 
B(x, r) is convex for every x € B. Therefore, for every x € B and for every pair of 
disjoint geodesic segments /3i,/?2 '• I '■= b) — > M contained in B(x,r), Theorem 
Eltells us that the function t i— > £(t) := d(f3i(t), fait)) is convex. If /3i,/?2 are not 
disjoint and neither of them is constant (in which case the result would be trivial) 
then we can only have the equality d((3\(t), fait)) = for a unique t = to, at which 
point the function £(t) attains an absolute minimum, and Theorem 1141 shows that 
£(t) is convex on (a, to) and on (to, b). But a real function which is convex on (a, to) 
and on (to, b), and which attains its minimum at to, must in fact be convex on all 
of 7= (a, b). 

This proves that the distance function d is convex on B(x, r) x B(x, r), for every 
x G B, which in turn means that d is uniformly locally convex near the diagonal. 

(2) In a Cartan-Hadamard manifold M every ball is convex, and two distinct 
geodesies in M can intersect in only one point (see ^| p. 259-261]), so the above 
argument applies globally. □ 

Remark 16. The assumption on curvature is necessary in order that d be uni- 
formly locally convex near the diagonal: it is easy to see that, for many disjoint 
nonconstant geodesic segments /3i and 02 in the sphere S 2 (take for instance two 
parallel meridians near the equator), the function t — > d(f3\(t), ^(t)) is not convex. 
Furthermore, as we will see in the next section, an important consequence of The- 
orem^] fails in the sphere S 2 , so the assumption on the jointly convexity of the 
distance function d : M x M — ► R near the diagonal seems to be much more than a 
mere technical requirement and is probably a necessary condition for the functions 
fx to be convex whenever / is. 

We conclude with a Corollary that sums up what the results we have shown tell 
us in the case of a Riemannian manifold of nonpositive curvature. 

When M is a complete finite-dimensional Riemannian manifold of nonpositive 
curvature we have the following result. Recall that a convex function / on a finite- 
dimensional Riemannian manifold M that only takes finite values is automatically 
continuous (see 0). 

Corollary 17. Let M be a complete finite- dimensional Riemannian manifold with 
sectional curvature K < 0. Let f : M — *■ R be a convex function. Then, for every 
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bounded open convex set U with compact closure U , there exists Xq > such that 
the functions fx : M — > R, defined by 

f x (x) = M{f(y) + l-d(x,y) 2 }, 

y£M ZA 

are convex and C 1 smooth on U for all A G (0, Ao). Moreover, 

(1) fx converges to f uniformly on U. 

(2) fx<f for all A > 0. 

(3) fx has the same inf and the same set of minimizers as f . 

(4) fx has the same symmetries as f ( that is, if f is invariant with respect to 
an isometry T : M — > M, then so is fx). 

Note that in this result we do not allow / to take infinite values. We are able to 
deal with functions / : M — ► R U {+00} when we furthermore assume that M is a 
Cartan-Hadamard manifold (that is, a simply connected Riemannian manifold of 
nonpositive curvature), either finite or infinite-dimensional. Also recall that Cartan- 
Hadamard manifolds enjoy the property that every two points can be connected by 
a minimizing geodesic (see |17| ). 

Corollary 18. Let M be a Cartan-Hadamard manifold (either finite-dimensional 
or infinite- dimensional). Let f : M — > RUj+oo} be a lower- semicontinuous convex 
function. Then the functions f\:M—* R, defined by 

fx(x)= M{f(y) + -Ld(x,y) 2 }, 
yeM ZA 

are convex and C 1 smooth on all of M for all A > 0. Moreover, 

(1) fx < f for all A > 0. 

(2) fx has the same inf and the same set of minimizers as f . 

(3) fx has the same symmetries as f does. 

(4) If f is uniformly continuous on bounded sets fx converges to f uniformly 
on bounded sets. 



5. Some applications 

If X is a Hilbert space (or more generally a reflexive Banach space), it is well 
known that for every closed convex subset C of X the distance function to C, that is, 
x 1 ► d(x, C) is convex and C 1 smooth away from C (even though C might not have 
a smooth boundary), and, as a consequence, every such C can be approximated by 
C 1 smooth convex bodies. We next show how the results proved above allow us 
to extend these two theorems to the class of Cartan-Hadamard manifolds (cither 
finite or infinite-dimensional), and we also note that this result completely fails in 
the sphere S 2 : there are closed convex sets C of arbitrarily small diameter in S 2 
such that x i— > d(x, C) is not convex on any neighborhood of C . 

Corollary 19. Let C be a closed convex subset of a Cartan-Hadamard manifold. 
Then the distance function to C , x 1— » d(x, C) = mi{d(x, y) : y G C} is C 1 smooth 
on M\C and, moreover, the function x i— > d{x, C) 2 is C 1 smooth and convex on 
all ofM. 
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Proof. Define / : M 

if x e C, 
otherwise . 

The function / is lower semicontinuous and convex on M. According to Theorem 
[TJ the function / A : M -> R, 

h{x) = inf{/(y) + ^d(x,y) 2 } = M{±d(x,y) 2 : y e C} = ^(x.C) 2 , 

is C 1 smooth and convex on M for all A > 0. By taking A = 1/2 we get that the 
squared distance function to C is C 1 smooth and convex on M . □ 

Definition 20. We say that a subset C of a Ricmannian manifold M is a C 1 
smooth convex body of M provided C is closed, convex, has nonempty interior, and 
dC is a one-codimensional C 1 smooth submanifold of M, 

Corollary 21. Let C be a closed convex subset of a Cartan-Hadamard manifold, 
and let U be an open subset of M with d(C,M \ U) > 0. Then there exists a C 1 
smooth convex body D of M such that C C D C U . 

Proof. Since d(C,dU) > we can take r — ^d(C,dU) and define D = {x G M : 
d(x,C) < r}. It is clear that C C D C J7, and D happens to be a C 1 smooth 
convex body because x i— > t2(ai, C) 2 is C 1 smooth and convex, and the derivative of 
the function C) 2 is nonzero at every point i£M\C(a convex function has a 
null derivative only at the points, if any, where it attains its minimum). □ 

The following example shows that Theorem ^] and the above Corollaries are 
false in general if we do not require that the manifold M has nonpositive sectional 
curvature. 

Example 22. Let M be the sphere x 2 + y? + z 2 = 1 in M 2 endowed with its usual 
Riemannian metric. Let C be a closed geodesic segment of diameter less than a 
number e with < e < 1. It is easy to see that the function d(-, C) 2 defined on M 

by 

d(x,C) 2 = inf d(x,y) 2 , 

is not convex on any open neighborhood of C, that is, Corollary El fails m M 
(hence so does Theorem II 20 . 

Finally, it should be noted that there is a strong link between the regularization 
method we have just presented and the following Hamilton- Jacobi partial differen- 
tial equation: 

' du(t,x) , 1 || du{t,x) || a _ n 
dt + 2 II dx Wx — u 

u{Q,x) = f{x), 

where u : [0, oo) x M — ► K, / : M — > R U {+oo}. If we assume that M is a finite- 
dimensional Cartan-Hadamard manifold and / is convex and lower-semicontinuous, 
then the function 

u(t, x) := inf {f(y) + ±-d(x, y) 2 } for t > 0, u(0, x) = f(x) 

y£M It 

is the unique viscosity solution of (*) (see [2] for the definition of viscosity solution 
to Hamilton- Jacobi equations on Riemannian manifolds). This is not very difficult 



(*) 
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to show directly. Alternatively, one can prove that Theorem 3.6 and Section 7.2 
of [H] remain true when K™ is replaced with a finite-dimensional Cartan-Hadamard 
manifold. 
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